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We present an analysis with improved sensitivity to the light charged Higgs 
{mH+ < mt — mb) searches in the top quark decays t — )• bH^ — )• b{T^h'r) + c.c. in the 
tt and single t/t production processes at the LHC. In the Minimal Supersymmetric 
Standard Model (MSSM), one anticipates the branching ratio B{H^ — )• t^Vt) — 1 
over almost the entire allowed tan f3 range. Noting that the r"*' arising from the 
decay — )■ r^z^r are predominantly right-polarized, as opposed to the r"*" from 
the dominant background — t- t~^i't, which are left-polarized, a number of 
H'^ /W~^ — )■ T'^Vr discriminators have been proposed and studied in the literature. 
We consider hadronic decays of the r^, concentrating on the dominant one-prong 
decay channel — )• p^Vr- The energy and pT of the charged prongs normalised to 
the corresponding quantities of the are convenient variables which serve as 
polariser. We use the distributions in these variables and several other kinematic 
quantities to train a boosted decision tree (BDT). Using the BDT classifier, and 
a variant of it called BDTD, which makes use of decorrelated variables, we have 
calculated the BDT(D)-response functions to estimate the signal efficiency vs. the 
rejection of the background. We argue that this chain of analysis has a high sensitiv- 
ity to light charged Higgs searches up to a mass of 150 GeV in the decays t — )• bH~^ 
(and charge conjugate) at the LHC. For the case of single top production, we also 
study the transverse mass of the system determined using Lagrange multipliers. 
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I. INTRODUCTION 

In many extensions of the standard model (SM), the Higgs sector of the SM is enlarged 
by adding an extra doublet of complex Higgs fields. After spontaneous symmetry breaking, 
one finds three neutral Higgs bosons {h, H, A) and a pair of charged Higgs bosons, 
These neutral and charged Higgs bosons have been searched for in high energy experiments, 
in particular, at LEP and the Tevatron. None of these Higgses have been seen so far, 
and upper limits exist on all of them [I]. We will concentrate here on the charged Higgs 
searches, in which the two key phenomenological parameters are the charged Higgs mass, 
rriHi, and tan /3, the ratio of the two vacuum expectation values, tan /3 = V2/V1. The searches 
for the are model-dependent, and the exclusion limits (expressed as a contour in the 
mH± - tan /3 plane) have to be taken together with the underlying model. For example, in 
the so-called two-Higgs-doub let-models (2HDM), a stringent limit exists on m/^± from the 
measured branching ratio for B — Xs'y and the NNLO estimates of the same in the SM, 
yielding m^i > 295(230) GeV at the 95% (99%) C.L., for almost the entire tan/3 values of 
interest [2]. This limit can be easily evaded in other models, in particular, in the minimal 
supersymmetric model (MSSM). 

Direct if^-searches are limited by the center-of-mass energy in e+e^ — > H^H^ annihi- 
lation processes, where they can be produced via s-channel exchange of a photon or a Z 
boson. These searches assume for the branching ratios B{H^ — )■ t^Vt) + B{H^ cs) = 1 
and hold for all values of E{H^ — t- r+z/,-). In the 2HDM framework, the cross section in the 
Born approximation depends only on mfj± (modulo the known couplings) and the present 
hmit is mH± > 79.3 GeV at 95% C.L. obtained at Ecm(e+e-) = 209 GeV from LEP [1]. 
The mass range mH+ < rrit — rrih has been searched in the process pp — tiX at the Teva- 
tron, followed by the decay t — > bH^ (and its charge conjugate). For example, Altonen 
et al. [3j have searched for the decay t — > bH^, followed by — )■ cs in 2.2 fb~^ of pp 
collisions at Ecmipp) = 1-96 TeV, obtaining upper limits on B{t — > bH^) between 0.08 and 
0.32 (95% C.L.), assuming B{H^ cs) = 1. In the MSSM, this probes only a very small 
tan/3 region, namely /3 < 1, which is not favoured by theoretical considerations [1]. The 
search for t — t- bH^, followed by — t- r+z/^ by Abazov et al. [5J in 0.9~^ of pp collisions 
at the Tevatron yield upper limits on B{t — )■ bH^) between 0.19 and 0.25 (95% C.L.) for 
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mH+ = 80 — 155 GeV and E{H^ — )• r^Ur) = 1. This excludes a small region (tan/3 > 35 
and mfj+ = 100 — 120 GeV) [6J. Thus, it is fair to conclude that the searches of the charged 
Higgses over a good part of the - tan /3 plane in the MSSM is a programme that still has 
to be carried out and this belongs to the LHC experiments. In anticipation, searches for the 
in pp collisions at E^ui = 7—14 TeV at the LHC have received a lot of attention [7HT^ . 
There are two regions, namely mH+ < rrit — rrib, which will be looked into in both the tt 
pair production and in single top (or anti-top) production in pp collisions, followed by the 
decays t — )■ bH'^ and — )■ r^z/^, and for mH± above the top quark mass, in which case 
production mainly takes place through the process gb — )■ tH~^, followed dominantly by the 
decay — tb. However, despite larger branching fraction, it may be hard to distinguish 
the — i- tb mode from the bckground. For large tan /3, the decay mode — t- t+z/t- 
becomes discernible. In this paper, we will concentrate on the light if^-scenario. 

The decay channel — + z/,- will play the key role in the searches of the light 
if^-bosons. The r"*" leptons arising from the decays — )■ r+z/^ and — )■ r^Ur are 
predominantly left- and right-polarised, respectively. Polarisation of the influences the 
energy distributions in the subsequent decays of the r^. Strategies to enhance the H^- 
induced effects in the decay t — )• b(W^, H^) — )■ bij^Ur), based on the polarisation of the r"*" 
have been discussed at length, starting from the pioneering work [T3I4T6] to the production 
and decays of a ti pair at the hadron colliders Tevatron and the LHC [I7H2Tj. Also the 
effects of the (QED and QCD) radiative corrections on such distributions in the dominant 
(one-charged prong) decay channels r+ — )■ tt+z/^, p+z/^, a^j^z/^ and l^v^Vj. have been worked 
out [22]. Following these studies, the construction of the r^-jet (as well as &-jet) are of central 
importance in if^-searches. We use the dominant single-charged-prong decay — )■ p^Vr 
as the polariser. As — )■ vr^vr'' is the dominant decay mode, the energy and transverse 
momentum of the vr^ in the r^-jet become quantities of main interest for our study. Likewise, 
the distribution in the angle ■?/', defined as 

COSIp = — 2 2 ' (IJ 

plays an important role in our analysis. Since the energy-momentum vectors of the 6-jet 
and the can be measured, this distribution is measurable at the LHC. We also note that 
this distribution is different from the conventional definition of the angle [23], in which 
the invariant mass m^^ is measured instead of m^^. The other distributions that enter in our 
analysis are listed in the next section. 

Having generated these distributions, characterising the signal t — t- bH^ — t- b{T^UT-) 
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b{p^h'r)i^T) and the background t — )■ bW~^ — )■ hij^Ur) — )■ b{p'^h'r)i^T) events, we use a tech- 
nique called the Boosted Decision Tree (BDT) - a classification model used widely in data 
mining [21] - to develop an identifier optimised for the t — > bH^ decays. In our calculation, 
we use both BDT and a variant of it called BDTD (here D stands for decorrelated) , where 
possible correlations in the input variables are removed by a proper rotation obtained from 
the decomposition of the square root of the covariance matrix, to discriminate the signal 
events from the large backgrounds. We recall that this technique has been successfully used 
to establish the single top quark production in pp collisions at the Tevatron [251 ES] (see [27] 
for details). Recently, we have applied this technique to a feasibility study of measuring 
the CKM matrix element \Vts\ from the decay t — )■ Ws at the LHC@14 TeV, and have 
estimated that a benchmark with 10% accuracy for this decay mode with a 10^ rejection of 
the background t — )■ Wb can be achieved with an integrated luminosity of 10 (fb)~^ [28] . 
We show in this paper that a similar BDTD-based analysis holds great promise in light-if^ 
searches at the LHC both in the pp — )■ tiX pair production and in the single top (or anti- 
top) production pp t/tX. Furthermore, we show that using a transverse mass definition, 
as suggested in [20], the process pp — )■ t/tX followed by the decays t bH~^ ,bW~^ , allows 
one to determine rather sharp Jacobian peaks for the transverse mass of the iJ^-bosons. 
The conventional definition of the transverse mass ^0], which was very helpful in the de- 
termination of the transverse mass of the bosons, is less suited for constructing the 
corresponding mass of the bosons. 

We note that an analysis using an iterative discriminat analysis method similar to the 
one presented here was carried out by Hesselbach et al. |;31j. In particular, detailed Monte 
Carlo comparisons of several variables incorporating the spin effects in charged Higgs boson 
production were presented to separate the tbH^ signal from the standard model tt back- 
ground both at the Tevatron {^/s = 1.96 TeV) and the LHC {^/s = 14 TeV). However, 
there are several significant differences in the two studies, such as the distribution in cos ijj 
(defined in eq. (1)), which plays an important role in our analysis. In addition, we have 
studied the case of single top production at the LHC, pp — )■ t/i + X, followed by the decays 
t — )■ b{H^ /W^ — > t+z/t-) +C.C., which was not considered in Ref. [3T] . 

This paper is organised as follows: In section 2, we analyse the process pp — t- tiX at the 
LHC, followed by the decay chains t — > bW^ , bH^, and the subsequent decays {H^, W^) — >• 
r+z/^, together with the BDTD-based analysis of the signal {t — )■ bH~^) and the SM decay 
background (t — )■ bW~^). The BDTD response functions are then used to work out the signal 
efficiency vs. the background rejection. In section 3, we repeat this analysis for the single 
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top (or anti-top) production pp — )■ t/tX at the LHC Section 4 contains a brief summary. 

II. tt PRODUCTION AND THE DECAY CHAINS t bW+/H+ h{T+Vr) AT 

THE LHC 

A. Production cross sections 

Theoretical predictions of the top quark production at the LHC have been obtained by 
including up to the next-to-next-to-leading order (NNLO) corrections in the strong coupling 
constant [32H3S] using modern parton distribution functions (PDFs) [HUl IH7] . Typical esti- 
mates for a{pp — )■ tlx) range from 874^33 pb for mt = 173 GeV and y/s = 14 TeV [38j to 
943±4(kinematics)l49(scale)±12(PDF) pb [34j- Compared to the ti production cross section 
at the Tevatron, this is larger by two orders of magnitude. The cross sections at the lower 
LHC energies, 7 and 10 TeV, have also been calculated jHU |3H], with a{pp — tiX) ^ 400 
pb at 10 TeV and about half that number at 7 TeV. Thus, for the top quark physics, the 
dividends in going from 7 to 14 TeV are higher by a good factor 4. 

B. Top quark decays t — )• b(W^ , H^) and charged Higgs decays — )• cs,T~^i'r 

Top-quark decays within the Standard Model are completely dominated by the mode 

t^b + W+ , (2) 

due to VJb = 1 to a very high accuracy. In beyond-the-SM theories with an extended Higgs 
sector, a light charged Higgs can also be produced via 

t^b + H+ . (3) 

The relevant part of the interaction Lagrangian is 



= ^yli^^*^^^ ^^'^ + 75) + S (1 - l,)}u, (p,)] 

+ ^^g^i/+K(p.)(l-75)«KP/)], (4) 

where A, B and C are model-dependent parameters which depend on the fermion masses 
and tan/3: 



A = mtCot(3, B = mbtain/3, C = m^tan/3. 



(5) 
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The decay widths of processes ^ and (|3| in the Born approximation are [T8] : 

(m^ + ml) + {m^^~ ml) ' - 2M^] , (6) 



647rM^mt 

pBorn _ A ^ 

X [[mf cot^ /3 + tan^ (3) [ml + ml - Mjj) - Amfml] , (7) 

where A (x, y, z) = + y"^ + — 2xy — 2xz — 2yz is the triangle function. The total top 
quark decay width in the Born approximation is obtained by adding the two partial widths 

ptot, Born pBorn , pBorn /oX 

QED corrections in the total decay width of the top quark are numerically small. The 0{as) 
QCD corrections were calculated in [321110] (see, also Ref. [H]) and have the form: 

ptot _ pBorn+QCD , pBorn+QCD /q\ 
J- t,RC — ^ t^bW + ^ t^bH ' V^) 

■pBorn+QCD -ptot.Born , r \ n '^^ f tr 4??^ \ 
^ t^b(W,H) - U + JW,h) , JW,H - ^ 1 ^ ^ j ■ 

Thus, in the branching ratio B{t — t- bH^), also this QCD correction drops out. However, 
radiative corrections coming from the supersymmetric sector to B{t — t- bH^) are rather 
important. They have been calculated in great detail in the literature, in particular for 
the MSSM scenario in [l^HH] . and can be effectively incorporated by replacing the 6-quark 



mass mb in the Lagrangian for the decay t — t- bH^ by the SUSY-corrected mass T/jCorrected _ 
mh/[l + Af,]. The correction is a function of the supersymmetric parameters and, for 
given MSSM scenarios, this can be calculated using the FeynHiggs programme |15]. In 
particular, for large values of tan/3 (say, tan /3 > 20)), the MSSM corrections increase the 
branching ratio for t bH^ significantly. This, for example, can be seen in a particular 
MSSM scenario in a recent update |16], from where we show B{t bH^) as a function of 
tan/3, calculated for mt = 175 GeV and various assumed values of the charged Higgs mass, 
indicated in Fig. [1} 

Since we are treating the case of the light charged Higgs, there are essentially only two 
decay modes which are important: — )■ r+z/^ and — )■ cs. The branching ratio of 
interest to us B{H^ — )■ t^i't) is given by [18j : 

J- H^rvr + J- H^cs 



H^cs 



327rM2 



(m^ cot^ (3 + ml tan^ /3) . 



w 
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FIG. 1: Branching ratio B{t — )• H^h) in MSSM as a function of tan/3 for the charged Higgs masses 
as indicated. (Figure taken from [46j). 



C. Event generation, trigger 

We consider in this section the process pp — )■ tiX, with both the t and t decaying into Wh. 
Our trigger is the leptonic decay W~ e^z/g or W~ — > h^p^. The other decays via 
— 7- t^Vt-. This makes up our main background. The signal events are generated in which 
one of the t or t decays via — )■ hH^ (or its charge conjugate W" — )■ bH^), see Fig. |2} 
The other t or t then decays leptonically, as in our trigger. In the Minimal Supersymmetric 
Standard Model (MSSM), for large tan/3 and mH+ < iTit, the branching ratio for the decay 
— )■ cs is small and one anticipates the branching ratio B{H^ — t- r+z/,-) 1. This is the 
parameter space in which the analysis reported here is valid. Noting that the r"^ arising from 
the decay — > r^z/^ are predominantly right-polarized, as opposed to the r"*" from the 
dominant background — ?■ r+z/^, which are left-polarized, a number of /W'^ — ?■ r+z/7- 
discriminators have been proposed and studied in the literature. We have used the dominant 
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FIG. 2: Feynman diagram for gg — )• tt, followed by the decay t — )• b{H^ — )• T^Vr) and t 
b{W- e-Ve). 

single-charged-prong decay p^i^r as the r+ polariser. Having set these branchings, we 

have generated 50K events for the process — )■ — > bW~^{bW~), with all of them decaying 
according to the chain described earlier, i.e., W~ — ?■ e^z^e and — ?■ r+z/^, with all the r's 
forced to decay into p + Vt (here and below, charge conjugates are implied). In calculating 
the required luminosity, we take into account the corresponding branching ratios, which are 
as follows [T] 

B{W+^e+Ve) = (10.75 ± 0.13)% , 
}3{W+^T+Ur) = (11.25 ± 0.20)% , 

B{H+ T+Ur) = 1.0 , 

i3(r+ p+Ur) = (25.5 ± 0.10)% . (11) 

We also generate the same number (50K) signal events, for each of the following charged 
Higgs masses: 171^+ = 90, 110, 130, 150 GeV. As for the background process, we force the r"*" 
to decay into p'^i'r 100% of the time. These events are generated using PYTHIA 6.4 [47j 
and for the decays of the r^, we use the programme called TAUOLA [48j to incorporate the 

polarization information on the decay distributions. 

We impose the following acceptance and trigger cuts: 

• \rie\ < 2.5, with i = e,T 

• \Vb,b\ < 2.5 
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• Pt, > 20 GeV 

• > 10 GeV 

• Pr. , > 20 GeV 

In order to discriminate the signal and background, we have studied a number of distri- 
butions, summarized below. 

• Distribution in the angle ip, defined in eq. [TJ This is defined for both the decay 
chains: t bW — h{TVr) b{pUr)i^T and t — > bH — > birVr) — ?■ b{piJr)'^T- Since the 
energy-momentum vectors of the 6-jet and the can be measured, this distribution 
is measurable at the LHC. We also note that this distribution is different from the 
conventional definition of the angle ip [23], in which the invariant mass m^^ is measured 
instead of m^f,. 

• Energy and pt of the 6-jets from the decays t — t- bW^ and t — > bH~^. 

• Energy and pr of the r"*" jets from the decays — > t^Uj. and — t- r+i^^, concen- 
trating on the single-charged-prong decays r+ — > p^Ur- 

• The ratio of the energy and px of the r"*" jets and their accompanying 6-jet. 

• As a measure of the r polarisation, we consider the fractional energy and transverse 
momentum of the single-charged prong (vr"*" in r"'"-jet). 

• For the case of single top production, we also study the transverse mass of the system 
determined using Lagrange multipliers |29j . 

• These distributions are used to train a boosted decision tree (BDT). Using the BDT 
classifier, and a variant of it called BDTD, which makes use of decorrelated variables, 
we have calculated the BDT(D)-response functions to estimate the signal efficiency 
vs. the rejection of the background. 

The strategy adopted by us to search for the decays t — )■ bH'^ is somewhat different 
from the traditional cut-based analysis, as, for example, reported in [8]. There the idea 
is to suppress the SM-background as much as possible, making use of additional variables, 
such as the missing Et, satisfying EJ^^^^ > 50 GeV. Our idea is, instead, to train a boosted 
decision tree classifier for both the signal and background events. Eventually, for a realistic 
analysis of the LHC data, we may have to reintroduce some of the cuts to suppress other 
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non-tt background, such as coming from the process pp — )■ + jets, which may also fake 
our signal. 



D. Details of the Analysis 

In Fig. [3] (right-hand frame), we show the cos^p distributions for the standard model (SM) 
process p + p ^ ti + X, followed by the decay chain t — )• bW h{rur) — )■ b{ph'r)i^T)- In the 
left-hand frame, we show the same distribution when one of the t 01 i decays via the chain 
t — 7- bH — )► b{TUr) — ^ b{ph'r)i^T) , for four different charged Higgs masses, as already stated 
in the previous section. For lower values of mH+, the cos?/^ distribution falls less steeply 
than the SM background. As mH+ increases, the cos ip distributions become steeper and are 
essentially confined to the negative values of cosip. This distribution then provides one of 
the discriminators to be fed to the BDTD analysis. 




FIG. 3: The distribution dN / d cos ip for the tt production as measured in the decay chain t — >• 
bW —7- b{TV-T) —7- b[pT>-T)v.j-) (right-hand frame), and in t — )• bH — t- b[TVr) — )• b{pi>r)i'T) for four 
different charged Higgs masses, as indicated on the figure (left-hand frame). 

In Fig. |4] (right-hand frames), we show the distributions in the energy of the &-jet, E{b)^ 
and the transverse momentum of the 6-jet, prib) from the SM process process pp — )■ ttX, 
followed by the decay chain discussed above. In the left-hand frames, the corresponding 
distributions are shown for the charged Higgs case. We remark that for the charged Higgs 
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case these distributions are softer than those from the SM due to the different hehcity struc- 
ture of the decays. This effect becomes stronger as 711^+ increases due to phase space. As a 
result, these distributions add to the discrimination power of the BDTD analysis. Note that 
these distributions reflect the event characteristics at the generation level. Obviously, due 
to the semileptonic decays of the 6-quark, and other detector effects, they will be modified. 
However, we expect that the dilutions due to these effects are sub-dominant. 
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FIG. 4: Distributions in the Energy of the 6-jet, E(b), and transverse momentum of the &-jet, 
Pxib) from the process pp — )• ttX, followed by the decay t — )• W~^b (right-hand frames), and the 
same distributions for the decay chain t — )• H^b with the four indicated charged Higgs masses 
(left-hand frame). 

In Fig. [5] (right-hand frames), we show the distributions in the energy of the r-jet, E{t — 
jet), and in the transverse momentum of the r-jet, ^^(r — jet) from the SM process, followed 
by the decay chain discussed above. In the left-hand frames, the corresponding distributions 
are shown for the charged Higgs case. In these distributions, the energy and pr-spectra of 
the T-jet coming from the charged Higgs decays are harder than those coming from the SM 
process, and this difference becomes more marked as mH+ increases. This complementary 
behaviour is expected for the same reason as discussed earlier for Fig. |4| again reflecting 
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the differences in helicity and phase space. It goes without saying that these distributions 
increase the discrimination power of the BDTD analysis. 




p^(T-jet) [GeV] P^(x-jet) [GeV] 



FIG. 5: Distributions in the Energy of the r-jet, E{t — jet), and transverse momentum of the 
T-jet, pt{t — jet) from the process pp — )• ttX, followed by the decay t — )• W^b (right-hand frames), 
and the same distributions for the decay chain t — >• H^b with the four indicated charged Higgs 
masses (left-hand frame). 

To make this effect more marked, we show the ratio of the energy and p-r-spectra involving 
the r- and 6-jets in Fig. [6j The SM distributions are shown in the right-hand frames, and 
those from the charged Higgs in the left-hand frames. These distributions show clearly the 
different shapes of the distributions SM vs. Higgs. For example, putting a lower cut on the 
ratios E{T—}et)/E{b) > 1 or pt{t —jet) / pT^b) > 1, most of the SM background is eliminated, 
whereas the charged Higgs-induced distributions surviving this cut are considerably larger, 
with the discrimination becoming stronger as mH+ increases. 

In Fig. [7} we show the distributions in the fractional energy of the single-charged prong 
(7r+ in r"'"-jet), E{'k)/E{t — jet), and in the transverse momentum of the single-charged 
prong, Pt(71")/pt(t— jet) from the SM process (right-hand frames) and those coming from the 
charged Higgs-induced process (left-hand frames) for mH+ = 90 GeV. As remarked earlier. 
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FIG. 6: Distributions in the ratio E{T—iet)/E(b) andpT(T— jet)/p7^(6) from the process — )• ttX, 
fohowed by the decay t — )• W^b (right-hand frames), and the same distributions for the decay chain 
t —7- H'^b with the four indicated charged Higgs masses (left-hand frame). 

we are using the dominant single-charged-prong decay r"*" — > p'^Ur as the polariser. As 
already noted in [H], the fractional energy distributions in z = Ea/Et, from the r-decay 
products r — A+ missing energy, the effect of the polarization is most marked for the 
decays r+ — > tt+z/,- and r+ — > p^Ur- This has been worked out in the coUinear limit, i.e., 
for E-r/iTLr ^ 1. Our variables differ from the one used in [r4J, in that we normalize to 
the visible r-energy and the visible Pt{t ~ jet), and not to the total r-energy. With our 
normalization, the vr "'"-energy measured in the decays r+ — t- tt+P^ will be a delta function, 
peaked at 1 in the variables shown in Fig. [7], and hence we concentrate on the decay chain 
r"*" — )■ p'^i't- These distributions also provide strong discriminants for the BDTD analysis. 

Briefly, the generated input is used for the purpose of training and testing the samples. 
We provide the input in terms of the variables discussed earlier for the signal {t — )■ hH^) 
and the background (t — )■ hW^), obtained with the help of a Monte Carlo generator. This 
information is used to develop the splitting criteria to determine the best partitions of the 
data into signal and background to build up a decision tree (DT). The separation algorithm 
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FIG. 7: Distributions in the fractional energy of the single-charged prong (vr"*" in r^-jet), 
E{tt)/E{t — jet), and in the transverse momentum of the single-charged prong, PTi'^)/PT{T — jet) 
from the pp — t- ttX, followed by the decay t — )• W~^b (right-hand frames), and the same distri- 
butions for the decay chain t — t- H^b with the four indicated charged Higgs masses (left-hand 
frame) . 

used in splitting the group of events in building up DT plays an important role in the 
performance. The software called the Toolkit for Multivariate Data Analysis in ROOT 
(TMVA) |19] is used for the BDT(D) responses in our analysis. The response functions for 
pp — )■ tiX at a center-of-mass-energy ^/s = 14 TeV at the LHC, followed by the background 
process t — > bW~^ (in shaded blue) and the signal t — t- bH^ (in shaded red) are shown 
in Fig. [Sj The four frames shown in this figure correspond to the charged Higgs masses 
mH+ = 90, 110, 130 and 150 GeV. As can be seen that the separation between the signal 
and the background increases as 771^+ increases. This improved separation as a function of 
mH+ will, however, be compensated to some extent by the decreasing branching ratio for 
the decay t — > bH~^ , as shown in Fig. [l] ^6] , obtained by using FeynHiggs [50] . 

The corresponding background rejection vs. signal efficiency curves from the process 
pp — )■ tiX calculated from the previous BDTD response at ^/s = 14 TeV are shown in Fig. M 
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FIG. 8: BDTD response functions for pp — )• ttX, with y/s = 14 TeV, foUowed by the decay 
t — >• W~^b (SM) and the decay chain t — )• H~^b, with the SM background (in shaded blue) and the 
charged Higgs signal process (in shaded red) for four different charged Higgs masses. 

for the four charged Higgs masses, as indicated on the frames. For a signal efficiency value 
of 90%, the background rejection varies between 50% and 90% as we move from 171^+ = 90 
GeV to mH+ = 150 GeV. 

In order to calculate the significance of our signal, we do the following simplified calcu- 
lation. We consider the less preferred case for tan/3 = 10, for which the branching ratio 
B{t H+b) in the MSSM shows a dip, with B{t H+b) ~ 0.02 for mH+ = 90 GeV (see. 
Fig. [T]). For the process pp — > tiX, the trigger is based on the decay t — )■ bW^ — )■ bi^ui, with 

= e"^,/i^, which has a summed branching ratio of about 0.2. Since, in the large-tan /3 
limit we are working, B{H^ — )■ r^Vr) — 1, and the r^-decay mode we are concentrat- 
ing on is r+ — )■ p'^Vri which has a branching ratio of 0.25, the product branching ratio 
t — )■ bH^ —7- b{T^Vr) b{p~^Ur)i^T = 5 X 10~^, which taking into account the trigger is 
reduced to 1.0 x 10^^. For an integrated luminosity of 10 (fb)~^ at y/s = 14 TeV, and 
summing over the charge conjugated modes yielding a factor 2, this yields 2 x 10^ signal 
events. For the background events, resulting from the production and the SM decays from 
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FIG. 9: SM background rejection vs. charged Higgs signal efficiency for the four charged Higgs 
masses indicated on the figure from the process pp — )• tiX, with y/s = 14 TeV. 

the process pp — ?■ ttX, the corresponding product branching ratio is 2.5%, which together 
with the trigger branching gives 5 x 10"'^, resulting in 10^ background events, where we have 
again taken into account the factor 2 from the sum of the charge conjugated states. Using 
the BDTD analysis, we get for a 50% signal efficiency, a background rejection of 90%. Thus, 
our estimated significance will be 

g ^ signal events ^ ^ j^qq ^^2) 

'\/ -^background events 

A more realistic calculation should consider a factor of 2 reduction due to the acceptance 
cuts, discussed in section A, as well as the efficiency to tag two 6-jets which is another factor 
of 2, and the efficiency of reconstructing a r — jet, estimated as 0.3 [7J. This amounts to a 
factor of about 10 reduction in both the number of signal and background events, resulting 
in a significance of about 30. This is high enough to take another factor 2 reduction due to 
various other cuts, which will be inevitable in a detector-based analysis taking into account 
non-tt backgrounds, not estimated here. Of course, this significance goes down as mH+ 
increases, keeping tan/3 fixed. Thus, for example, for tan/3 = 10 and mH+ = 150 GeV, 
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the reduction in the number of events will be approximately 5 (a factor 10 decrease in 
— 7- H^b), compensated by a factor 2 increase in the signal efficiency calculated from 
the BDTD response). This would yield S 6, which is just above the discovery limit for a 
charged Higgs below the top quark mass. 

A number of checks has been performed in order to test the robustness of the results. For 
instance, the cut on the minimum transverse momentum of the r-jet has been raised from 
10 GeV to 20 GeV. The corresponding figure displaying the background rejection vs. the 



charged Higgs signal efficiency is shown in Fig. 10 A comparison with Fig. |8| obtained 
with a 10 GeV cut on the minimum transverse momentum of the r-jet, shows that the two 
figures are very similar. The price to pay for the acceptance is, relatively speaking, minor, 
going down from 0.6 to 0.5. We had conservatively taken this to be 0.5 in our numerical 
calculations. 
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FIG. 10: SM background rejection vs. charged Higgs signal efficiency for the four charged Higgs 
masses indicated on the figure from the process pp — )• tiX. This figure differs from the one shown 
in Fig. [9] in the minimum transverse momentum of the r-jet, which is set to 20 GeV as opposed to 
10 GeV used in the earlier figure. 

The above analysis presented for the LHC energy ^/s = 14 TeV has been repeated for a 
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center of mass energy a/s = 7 TeV, at which energy the LHC is collecting data currently. As 
of preparing this report, the integrated luminosity of the LHC is above 1 inverse femtobarn, 
and the projection for end 2012 is of order 10 inverse femtobarns. We have generated events 
at a/s = 7 TeV, and have calculated all the distributions presented earlier for 14 TeV. The 
shapes of these distributions are essentially similar. This is reflected in the BDTD response 



functions for the SM background and the charged Higgs signal, presented in Fig. 11, and 



in the SM background rejection vs. the charged Higgs signal efficiency, shown in Fig. 12 



However, the cross sections for pp tiX at 7 TeV is approximately a factor 4 smaller 
than at 14 TeV [SU |38]. This implies that our calculations for the significance obtained 
at ^/s = 14 TeV have to be divided by a factor 2 to get the corresponding significance at 
^/s = 7 TeV. This will reduce the sensitivity of the charged Higgs in tan /3- mH+ plane. 
For example, for mH+ close to the kinematic limit rrit — mb, a signal is expected only for 
tan/3 > 20. 
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FIG. 11: BDTD response functions for pp — t- ttX , with ^/s = 7 TeV, followed by the decay 
t — >• W~^b (SM) and the decay chain t — >• H~^b, with the SM background (in shaded blue) and the 
charged Higgs signal process (in shaded red) for four different charged Higgs masses. 

A potential dilution of the polarization information has to be kept in mind. The single 



19 





FIG. 12: SM background rejection vs. charged Higgs signal efficiency for tlie four charged Higgs 
masses indicated on the figure from the process pp — )• tiX with ^/s = 7 TeV. 

charged-prong hadronic decays of the are essentially made up of the decays — )■ vr^ + Ur 
(with a branching ratio of 10.9%), — > > vr^vr^) + Ur decays (with a branching ratio of 
25.5%) and — )■ cbf{—^ tt^tt^tt^) + u^- decays (with a branching ratio of 9.3%). Separating 
the TT^Ur mode from the p^i^r niode should, in principle, be possible due to the lack of 
deposited energy in the 7r° or electromagnetic cluster accompanying the vr^ in the former, 
but separating the p^z/,- mode from the afu^- mode will not be easy. Fortunately, the 
branching ratio of the latter is only 40% of the former. So, the number of 7r° clusters (0, 1 
and 2) will have to be included in the analysis as a new variable. The — )■ p^Uj. decays 
dominate the one charged track {ir^jK^) and one electromagnetic or 7r° cluster. However, 
we stress that the BDT can be trained to reduce the dilution. 

In a realistic analysis, a further source of reduction in our estimates of the significance 
would come from the wrong assignment of the 6-jet charges, though this effect is minor 
compared to the ones discussed above. The &-jet charge identification efficiency is estimated 
at present to be around 65% [7] , using standard techniques based on a weighted average of the 
charges of the particles in the jet, with the weights being proportional to their momenta. 



20 



However, a simple algorithm can be designed, which takes into account in addition the 
angular correlations between the trigger lepton, the tau-jet and the charges, reducing the 
6-jet mis-assignment to about 20% for the charged Higgs masses close to the mass. For 
higher charged Higgs masses, this can be further brought down by simply taking the 6-jet 
with the smaller (larger) transverse momentum to be that associated with the charged Higgs 
(resp. W^) boson. 

We also mention that we have not considered the background from the process pp ^ ti ^ 
{bii'i){bjj). However, it has been shown in [7j that this background can be well separated 
in a standard cut analysis from the — )■ — )■ {Mvi){hTVr) process. With our TMVA 
approach, this background will be tamed though we will have to introduce also the missing 
Et as a variable in the BDT training. We are aware of the non-tt background, which are 
dominated by the Z + jets and W + jets. These have been studied in great detail in [8], 
with the conclusion that they can be brought below the signal by the additional use of the 
E™^^^-cni. We have not used the E^^^^-cvX, as we have concentrated only on the SM tiX 
background, but will do so in a more realistic detector-based analysis in the future. 



III. SINGLE t/t PRODUCTION AND THE DECAY CHAINS 

t bW+/H+ b{T+Ur) AT THE LHC 

A. Cross sections at the LHC 



The single top (or anti-top) cross sections in hadron hadron collisions have been calculated 
in the NLO approximation [HIHSSj. Recalling that there are three basic processes at the 
leading order which contribute to a{pp — )■ t/iX), namely the t-channel: qb — )■ q't, the s- 
channel: qq' — > bt; and the associated tW production bg — )■ tW~ , the cross section estimated 
at the Tevatron is a{pp — > tX) = cr{pp — )■ iX) ^ 1.8 pb for both the top and anti-top 
production. At the LHC@14 TeV, one estimates a{pp — ?■ tX) ~ 200 pb and about half this 
number for a{pp — > iX), yielding the summed single top and anti-top cross sections at about 
300 pb, also approximately two orders of magnitude larger than those at the Tevatron. With 
a luminosity of 10 fb~^, one anticipates 0(3 x 10^) single top (or anti-top) events. 

As mentioned in the introduction, there are three different mechanisms of producing 
a single top (or anti-top) quark in hadronic collisions, the s-channel, the t-channel, and 
the associated production tW-channel. The Feynman diagram for the dominant t-channel 



partonic process qb — > q't, followed by the decay t — > b{H^ — > r+i/^) is shown in Fig. 13 
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The partonic cross section is then convoluted with the parton distribution functions to 
calculate the cross sections in pp t + X and pp ^ t + X. Since, we are using PYTHIA 
6.4 |17j to do the simulation of the single top (or anti-top) production, not all channels 
are encoded there yet. However, as we use the generator to calculate the acceptance only, 
but the total cross sections are normalized to the theoretical calculations, the estimates 
presented here should hold approximately. Since most of the distributions calculated by us 
for the processes pp — )■ tix and pp — )■ t/iX are in the same variables, we comment only 
briefly on the distributions for the signal t — )■ bH^ — hr^Vr and the background process 
t bW+ br+Ur. 




FIG. 13: Feynman diagram for qb — )• q't, followed by the decay t — )• b{H^ — )• r^Ur). 



In Fig. 14, we show the distribution dN/dcosip for the pp — )■ t/i + X production as 
measured in the decay chain for the SM background process t — )■ bW — )■ bijUT-) b{pUr)^T) 
(right-hand frame), and for the signal t — t- bH — )■ b{TVT-) — )■ b{pi'r)^T) (left-hand frame) for 
four different charged Higgs masses, as indicated on the figure. The SM background in the 
process pp — )■ t/i + X falls more steeply as a function of cosip than is the case for the ti 
production pp ^ ti + X, due to the acceptance cuts. The trend is similar in the signal 
process. However, also in the single top (or anti-top) production, this distribution provides 
a good discriminant as input to the BDTD analysis. 

The distributions in the energy of the 6-jet, E{b), and transverse momentum of the b- 
jet, prib) from the process pp — t/iX, followed by the SM decay t — )■ W^b are shown in 



Fig. 15 (right-hand frames), and the same distributions for the decay chain t — )■ H^b with 
the four indicated charged Higgs masses are also shown in this figure (left-hand frame). 
These distribution are very similar to the ones shown for the pp — ?■ tiX processes, as they 
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FIG. 14: The distribution dN /d cos ip for the pp ^ t/t + X production as measured in the decay 
chain t — )• hW — )• b{Ti>r) — )■ b{pi'r)i^T) (right-hand frame), and in t — )• bH — )• b{Ti'r) — ^ b{pi>T)vr) for 
four different charged Higgs masses, as indicated on the figure (left-hand frame). 

essentially reflect the kinematics of the decays t — )■ W'^h and t — )■ H^h. 



In Fig. 16, we show the corresponding distributions for the r-jet, E{t — jet), and for 
the transverse momentum of the r-jet, Pt{t ~ jet) from the process pp — > t/iX, followed 
by the SM decay t — > W^h (right-hand frames), and the same distributions for the decay 
chain t — > H^h with the four indicated charged Higgs masses (left-hand frame). These 
distributions, likewise, are very similar to the ones shown for the ti production case, shown 
in the previous section. 

The distributions in the ratio E{t — ]ei) / E{h) and Pt(t ~ jet)/pT(&) from the process 



pp — t/tX, followed by the SM decay t — )■ W^h are shown in Fig. 17 (right-hand frames), 
and the same distributions for the decay chain t — t- H^h are also shown in this figure with the 
four indicated charged Higgs masses (left-hand frame). As anticipated, these distributions 
are also similar in the single top (anti-top) production and in the ti production. 

The effects of different chiralities in the SM decay chain t — hW^ — t- hij^Vr) followed 
by the r+ decay r+ — )■ p^Ur, and in the signal process t — )■ bH^ 6(r+z/T-) followed by the 
r+ decay r+ — t- p^Ur are shown in Fig. [Tsj Once again, these distributions in the fractional 
energy of the single-charged prong (tt"*" in r+-jet), E{ti)/E{t — jet), and in the transverse 
momentum of the single-charged prong, Pxi^^) / Pt{'T' ~ jet) are very similar in the processes 
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FIG. 15: Distributions in the energy of the b-jet, E{b), and transverse momentum of the 6-jet, 
Pxib) from the process pp — )■ t/iX, followed by the decay t — >• W~^b (right-hand frames), and the 
same distributions for the decay chain t — )■ H^b with the four indicated charged Higgs masses 
(left-hand frame). 

PP ~^ t/iX and pp — > tiX, as expected. 

One important difference between the analysis of the single top (or anti-top) production 
compared to the tt production process lies in the fact that the missing transverse energy 
and momentum can be ascribed in the former to the r-neutrino, Ur- This is different in the 
case of the ti production, as one of the t-oT t-quarks decays via t — )■ bW~^ — )■ bi^ue, which is 
used as a trigger. Thus, the missing transverse energy or momentum can not be traced to 
the decay of the r-lepton alone in the case of ti production. As already stated in [29| the 
missing transverse energy and momentum profile in the case of the single top (or anti-top) 
process pp — )■ t/tX followed by t — )■ bH^ — )■ br^Vj. can be used to constrain the mass of 
the charged Higgs. We pursue this idea, by using two different definitions of the transverse 
mass. In the first case, called m'j} , this is defined as in |30j : 

= 2p^p^(l — cos (ply) , (13) 
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FIG. 16: Distributions in the energy of the r-jet, E{t — jet), and transverse momentum of the 
r-jet, pt{t — jet) from the process pp — t/tX, foUowed by the SM decay t — W^b (right-hand 
frames), and the same distributions for the decay chain t — t- H^b with the four indicated charged 
Higgs masses (left-hand frame). 

where p^, p^, and (pii, are the momenta and angle between the leptons in the plane perpen- 
dicular to the pp collision axis. This definition was proposed to determine the transverse 
mass of the boson in pp collisions using the decay modes — )• e^z/g and — )• /x^z/^. 
In our case, where the charged Higgs decays via — )■ r the charged lepton is the r"*", 
which is not measured experimentally. Since, we use the decay r+ — )• p^Ur, we replace the 
p^ by the px of the p+. The resulting ■'-distributions are shown in the upper two frames 



in Fig. 19 for the SM background (right-hand frame) and the charged Higgs case (left-hand 
frame). As seen from the distributions shown in the left-hand frame, this definition is not 
useful to see the Jacobian peak in the transverse mass of the H^. This is anticipated since 
there are are two undetected neutrinos from the vertex. The distributions in m^"* for the 
SM (IV^)-background and the if^-signal are different, and they do add to the discriminating 
power in the BDTD analysis. 

For the processes pp{p) — ttX and pp{p) — > tbX, with the subsequent decay of t — )■ 
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FIG. 17: Distributions in the ratio E{t — jet)/ E{b) and ~ i^^)/PT(b) from the process 

PP ~^ t/tX, fohowed by the SM decay t — >• W~^b (right-hand frames), and the same distributions 
for the decay chain t — t- H^b with the four indicated charged Higgs masses (left-hand frame). 

b{W^ , H^), if one of the two b jets could be associated with the semileptonic decay of 
the top quark, then the on-shell constraint for the top quark could be used in the form 
^^miss + Pb)'^ = Tn"^. In this case, a transverse Higgs mass can be defined by maximizing 
the invariant mass, {m!j^Y = max[(p£ + p^^^^)'^], since it is bounded from above by the top 
quark mass, with the charged Higgs transverse mass satisfying 'mH+ < friff < rrit, where 
mH+ is the true chraged Higgs mass. This leads to the following transverse mass definition 
for [29], which we call 

« ? = (v/^^ + IpI+pI+pT")' - Pt) - ifr + P^n' ■ (14) 

(2) 

This expression holds by neglecting the 6-quark mass. We have calculated the distribu- 
tions, by replacing the p^ (which is p^ for our case) by p^. These distributions are shown in 



the lower two frames of Fig. 19 with the SM background (yielding the Jacobian peak of the 
W^) shown on the right-hand frame, and the corresponding Jacobian peaks for the charged 
Higgs case, shown in the left-hand frame. For all the four charged Higgs masses shown in 
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FIG. 18: Distributions in the fractional energy of the single-charged prong (vr''" in r^-jet), 
E{tt)/E{t — jet), and in the transverse momentum of the single-charged prong, PTi'^)/PT{T — jet) 
from the pp — t- t/tX, followed by the SM decay t — )• W'^b (right-hand frames), and the same 
distributions for the decay chain t — t- H^b with the four indicated charged Higgs masses (left-hand 
frame) . 

this frame, the Jacobian in m)p has a sharp peak. Measuring these distributions provides, 
in principle, an estimate of H^. We will use these distributions in to train our BDTD 
sample. 

The distributions generated and discussed have been used to train the BDTD algorithms 



and the resulting response functions are shown in Fig. 20 The separation between the 
signal and the background improves as 171^+ increases, a trend which was also observed in 
the pp — 7- tiX production process. 

The corresponding background rejection vs. signal efficiency curves from the processes 
PP ~^ t/iX calculated from the previous BDTD response at ^/s = 14 TeV are shown in 
Fig. [21] for the four charged Higgs masses, as indicated on the frames. For a signal efficiency 
value of 90%, the background rejection varies between 40% and 99% as we move from 
mH+ = 90 GeV to mH+ = 150 GeV. 
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FIG. 19: Transverse mass distributions for the in the process pp — t- t/i + X followed by the 
decay t — t- bW~^ (right-hand frames) and for the transverse mass for the decay chain t — t- H~^b 
with the four indicated charged Higgs masses (left-hand frame). The definitions used for defining 
the transverse masses and are given in the text. 

In order to calculate the significance of our signal, we do the following simplified calcula- 
tion. We consider again the less preferred case for tan /3 = 10, with B{t — > H^b) ~ 0.02 for 
mH+ = 90 GeV. For the process pp — > t/tX, our trigger is based on the r-jet, coming from the 
decays W^/H^ — > r^z/,-. Since, in the large-tan (3 limit we are working, B{H^ — t- t^Ut) — 1, 
and the r^-decay mode we are concentrating on is r+ — t- p^Ur, which has a branching ratio 
of 0,25, the product branching ratio t — t- bH^ — t- 6(r+z/T-) — )■ b{p^Ur)i^T = 5 x 10^'^, which is 
the same as in the case of top-quark pair production process pp — )■ tiX. For an integrated 
luminosity of 10 (fb)~^, and inclusive single top cross section a{pp t/tX = 300 pb at 
^/s = 14 TeV, this yields 1.5 x 10^ signal events. For the background events, resulting from 
the production and the SM decays from the process pp — ?■ t/tX, the corresponding product 
branching ratio is 2.5%, again the same as in the case of top-quark pair production process 
pp — )■ tiX. yielding 7.5 x 10^ background events. Using the BDTD analysis, we get for a 
50% signal efficiency, a background rejection of 90%. Thus, our estimated significance will 



28 



' I ' ' 

]mj^ = 90 GeV 
aSM Background 



-0.4 




0.2 0.4 0.6 
BDTD 



im„ = 130 GeV 
I SM Background 




"T" 



0.2 0.4 0.6 0.8 
BDTD 




0.2 0.4 0.6 
BDTD 




FIG. 20: BDTD response functions for pp — t- t/tX, followed by the decay t — )■ W^b (SM) and 
the decay chain t — )• H~^b, with the SM background (in shaded blue) and the charged Higgs signal 
process (in shaded red) for four different charged Higgs masses. 



be 



S 



signal events 



7.5 X 10^ 



- -1 85 . (15) 

Y -^background events \/7.5 X 10"^ 

A more realistic calculation should consider a factor of 2 reduction due to the acceptance 
cuts, discussed in section A, as well as the efficiency to tag the &-jet, estimated as 70%, 
and the efficiency of reconstructing a r — jet, estimated as 0.3. This amounts to a factor 
of about 10 reduction in both the number of signal and background events, resulting in a 
significance of about 25. Of course, this significance goes down as mii+ increases, keeping 
tan/3 fixed. Thus, for example, for tan/3 = 10 and mf{+ = 150 GeV, the reduction in the 
number of events will be approximately 5 (a factor 10 decrease in B{t — ?■ H^b), compensated 
by a factor 2 increase in the signal efficiency calculated from the BDTD response). Since 
the background rejection goes up to 99%, this would yield S ~ 25, allowing to search for a 
charged Higgs in the decay t — )■ bH^, essentially up to a charged Higgs mass close to the 
kinematic limit. 
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FIG. 21: SM background rejection vs. charged Higgs signal efficiency for tlie four charged Higgs 
masses indicated on the figure from the process pp — )■ t/tX. 

We would like to stress that our philosophy in this paper is to show how to disentangle 
the process — )■ t + X — > H'^b + X from pp ^ t + X ^ W^h + X. In particular, single 
top production in hadron colliders is subject itself to backgrounds [7] which we have not 
considsered here. The most relevant of these backgrounds is the Whh production. Needless 
to say that the cosip, the polarisation information on the from the decay — )■ p^Ur, 
and the transverse mass distribution will retain their discriminant power to suppress them, 
albeit at the cost of a small loss in the significance of the signal. We plan to take this into 
account together with a complete treatment of the detector effects in a forthcoming more 
realistic analysis, which is required to assign an error on the charged Higgs mass due to such 
effects. 



IV. SUMMARY AND OUTLOOK 

We have reported here an analysis with improved sensitivity to charged Higgs searches in 
top quark decays t — > bH^ — )■ hr^Vr at the LHC We concentrate on hadronic decays, in 
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particular, the decay mode r — ?• p z/,-, and take into account the polarisation information 
of the passed on to . The observables which play a dominant role in our analysis are 
the energy and pt of the 6-jets from the decays t — > hW^ and t — )■ feif"*", energy and px of 
the r^-jets from the two decay chains, and the energy and pt of the single-charged prong 
(tt^ coming from the decay chain — )■ p^z/^ — )■ tt^z/,-). Distributions in these variables are 
studied together with angular distribution in cos ip defined in eq. [T} This information is fed 
to a multivariate analysis using the BDTD techniques. The BDTD response shows that a 
clear separation between the t — hW'^ and t hH^ decays can be achieved in both the 
tiX pair production and the t/iX single top production at the LHC. We have also shown 
that using a transverse mass definition, as suggested in [22], the process pp — )■ t/tX allows 
one to determine sharp Jacobian peaks for the mass of the if^-bosons. With the benchmark 
integrated luminosity of 10 fb~^ at 14 TeV, the light charged Higgs (m//+ < rrit — rrib) can 
be discovered for all values of tan/3, where the decay mode — )■ r^Uj. is dominant. 

In estimating the quoted significances, we have assumed that the decay t — > bW~^ 
makes up the dominant background. This should be refined by taking into account non-t- 
backgrounds, such as coming from {Z,W) + jets. 

Acknowledgements: We thank Merlin Kole, Theodota Lagouri and Torbjorn Sjostrand 
for helpful discussions. This research was partially supported by MICINN (Spain) under 
contract FPA2008-00601. 



[1] K. Nakamura et al. [ Particle Data Group Collaboration ], J. Phys. G G37, 075021 (2010). 



[2] M. Misiak et al, Phys. Rev. Lett. 98, 022002 (2007) |arXiv: hep-ph/0609232|. 
[3] T. Aaltonen et al. [CDF Collaboration], Phys. Rev. Lett. 103, 101803 (2009) [arXiv:0907.1269l 
[hep-ex]]. 

[4] For a recent update, see G. Bernardi, M. Carena, and T. Junk, "Higgs Bosons: Theory and 

Searches", in the PDG review [1]. 
[5] V. M. Abazov et al. [DO Collaboration], Phys. Rev. D 80, 051107 (2009) jarXiv: 0906. 5326] 

[hep-ex]]. 



[6] O. Stal, arXiv:1012 .2709 [hep-ph]. 

[7] G. Aad et al. [The ATLAS Conaboration], |arXi v:0901.0512| [hep-ex]. 

[8] ATLAS Cohaboration, ATLAS-NOTE ATL-PHYS-PUB-2010-006 (2010). 

[9] ATLAS Collaboration, ATLAS-NOTE ATL-PHYS-PUB-2010-003 (2010). 



31 



[10] G. L. Bayatian et al. [CMS Collaboration], J. Phys. G 34 (2007) 995. 

[11] M. Baarmand, M. Hashemi and A. Nikitenko, J. Phys. G 32 (2006) N21. 

[12] R. Kinnunen, PoS CHARGED2008, 007 (2008). 

[13] K. Hagiwara, A. D. Martin and D. Zeppenfeld, Phys. Lett. B 235, 198 (1990). 

[14] B. K. Bullock, K. Hagiwara and A. D. Martin, Phys. Rev. Lett. 67, 3055 (1991). 

[15] A. Rouge, Z. Phys. C 48, 75 (1990). 

[16] B. K. Bullock, K. Hagiwara and A. D. Martin, Nucl. Phys. B 395, 499 (1993). 

[17] D. P. Roy, Phys. Lett. B 277, 183 (1992). 

[18] S. Raychaudhuri and D. P. Roy, Phys. Rev. D 52, 1556 (1995) |arXiv:hep-ph/ 9503251]; Phys. 



Rev. D 53, 4902 (1996) l| arXiv:hep-ph/9507388| . 
[19] D. P. Roy, Phys. Lett. B 459, 607 (1999) |arXiv:hep-ph/9905542| . 
[20] K. A. Assamagan and Y. Coadou, Acta Phys. Polon. B 33 (2002) 707. 

[21] Q. H. Cao, S. Kanemura and C. P. Yuan, Phys. Rev. D 69, 075008 (2004) |arXiv:hep- 
|ph/03li083i . 



[22] A. Ah, E. A. Kuraev and Y. M. Bystritskiy, Eur. Phys. J. C 67, 377 (2010) | ;arXiv:0911.3027 
[hep-ph]]. 



[23] D. Eriksson, G. Ingelman, J. Rathsman and O. Stal, JHEP 0801, 024 (2008) |arXiv:0710.5906 
[hep-ph]]. 

[24] J. Han and M. Kamber, Data mining concepts and techniques, Elsevier 2006. 



[25] V. M. Abazov et al. [DO Collaboration], Phys. Rev. Lett. 103, 092001 (2009) [arXiv:0903.6850 
[hep-ex]]. 



[26] T. Aaltonen et al. [CDF Cohaboration] , Phys. Rev. Lett. 103, 092002 (2009) [ arXiv:0903.0885 
[hep-ex]]. 

[27] Z. Liu, FERMILAB-THESIS-2009-45. 

[28] A. Ah, F. Barreiro, T. Lagouri, Phys. Lett. B693, 44-51 (2010). |arXiv:1005.4647| [hep-ph]]. 



[29] E. Gross, O. Vitells, Phys. Rev. D81, 055010 (2010). |arXivr0907.5367 [hep-ph]] 



[30] J. Smith, W. L. van Neerven and J. A. M. Vermaseren, Phys. Rev. Lett. 50, 1738 (1983). 
[31] S. Hesselbach, S. Moretti, J. Rathsman and A. Sopczak, Eur. Phys. J. C 53, 311 (2008) 



,arXiv:0708. 4394 [hep-ph]]. 
[32] R. Bonciani, S. Catani, M. L. Mangano and P. Nason, Nucl. Phys. B 529, 424 (1998) [Erratum- 

ibid. B 803, 234 (2008)] |arXiv:hep-ph/9801375l . 
[33] M. Cacciari, S. Frixione, M. L. Mangano, P. Nason and G. Ridolfi, JHEP 0809, 127 (2008) 

[arXiv:0804.2800 [hep-ph]]. 



32 



[34] N. Kidonakis and R. Vogt, Phys. Rev. D 78, 074005 (2008) [ arXiv: 0805 .3 844' [hep-ph]]. 

[35] S. Moch and P. Uwer, Phys. Rev. D 78, 034003 (2008) |a rXiv:0804.l"476 | [hep-ph]]. 

[36] A. D. Martin, W. J. Stirhng, R. S. Thorne and G. Watt, Phys. Lett. B 652, 292 (2007) 

|arXiv:0706.0459l [hep-ph]]. 
[37] P. M. Nadolsky et al, Phys. Rev. D 78, 013004 (2008) [arXiv:0802.0007| [hep-ph]]. 



[38] U. Langenfeld, S. Moch and P. Uwer, arXiv:0907.2527 [hep-ph] 



[39] A. Czarnecki and S. Davidson, Phys. Rev. D 47, 3063 (1993) [arXiv:hep-ph/9208240] . 



[40] A. Czarnecki and S. Davidson, Phys. Rev. D 48, 4183 (1993) [arXiv:hep-ph/9301237] . 
[41] C. S. Li and T. C. Yuan, Phys. Rev. D 42, 3088 (1990) [Erratum-ibid. D 47, 2156 (1993)]. 



[42] J. Guasch, R. A. Jimenez and J. Sola, Phys. Lett. B 360, 47 (1995) | arXiv:hep-ph/9507461] . 
[43] J. A. Coarasa Perez, D. Garcia, J. Guasch, R. A. Jimenez and J. Sola, Eur. Phys. J. C 2, 373 

(1998) | arXiv:hep-ph/9607485 ] . 
[44] M. S. Carena, D. Garcia, U. Nierste and C. E. M. Wagner, Nucl. Phys. B 577, 88 (2000) 



|arXiv:hep-ph/9912516| . 

[45] T. Hahn, S. Heinemeyer, W. Hollik, H. Rzehak and G. Weiglein, Nucl. Phys. Proc. Suppl. 
183, 202 (2008). 



[46] A. Sopczak, PoS CHARGED2008, 023 (2008). |arXiv: 0907.1498 [hep-ph]]. 

[47] T. Sjostrand, S. Mrenna, P. Z. Skands, JHEP 0605, 026 (2006). [hep-ph/OeOSTTS] . 

[48] S. Jadach, Z. Was, R. Decker and J. H. Kiihn, Comput. Phys. Commun. 76, 361 (1993). 

[49] A. Hocker et al, PoS A CAT, 040 (2007) [ar Xiv:phy sics/0703039|. 

[50] S. Heinemeyer, W. Holhk, G. Weiglein, Comput. Phys. Commun. 124, 76-89 (2000). | hep-| 
|ph/9812"320l . 

[51] B. W. Harris, E. Laenen, L. Phaf, Z. Sullivan and S. Weinzierl, Phys. Rev. D 66, 054024 



(2002) tarXiv:hep-ph/0207055] . 



[52] N. Kidonakis, Phys. Rev. D 74, 114012 (2006) |arXiv:hep-ph/0609287]. 



[53] Q. H. Cao, R. Schwienhorst and C. P. Yuan, Phys. Rev. D 71, 054023 (2005) |arXiv:hep- 
[ph/0409040|. 

[54] S. Heim, Q. H. Cao, R. Schwienhorst and C. P. Yuan, Phys. Rev. D 81, 034005 (2010) 



arXiv:091 1.0620 [hep-ph]]. 



[55] N. Kidonakis, Phys. Rev. D 75, 071501 (2007) [ arXiv:hep-ph/0701080| . 
[56] N. Kidonakis, Nucl. Phys. A 827, 448C (2009) jarXiv: 090 1.21551 [hep-ph]]. 



